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ABSTRACT 

From optical R band data of the double quasar QSO 0957+56lA,B, we made two 
new difference light curves (about 330 days of overlap between the time-shifted light 
curve for the A image and the magnitude-shifted light curve for the B image). We 
observed noisy behaviours around the zero line and no short-timescale events (with a 
duration of months), where the term event refers to a prominent feature that may be 
due to microlensing or another source of variability. Only one event lasting two weeks 
and rising - 33 mmag was found. Measured constraints on the possible microlensing 
variability can be used to obtain information on the granularity of the dark matter in 
the main lensing galaxy and the size of the source. In addition, one can also test the 
ability of the observational noise to cause the rms averages and the local features of the 
difference signals. We focused on this last issue. The combined photometries were re- 
lated to a process consisting of an intrinsic signal plus a Gaussian observational noise. 
The intrinsic signal has been assumed to be either a smooth function (polynomial) or 
a smooth function plus a stationary noise process or a correlated stationary process. 
Using these three pictures without microlensing, we derived some models totally con- 
sistent with the observations. We finally discussed the sensitivity of our telescope (at 
Teide Observatory) to several classes of microlensing variability. 

Key words: gravitational lensing - dark matter - galaxies: elliptical and lenticular, 
cD - quasars: individual: Q0957+561 



1 INTRODUCTION 

1.1 Microlensing caused by the Galaxy and other 
spirals 

Dark matter dominates the outer mass of the Milky Way. 
In principle, the population of the Galactic dark halo may 
include astrophysical objects as black holes, brown dwarfs, 
cool white dwarfs, etc., i.e., MACHOs (massive compact halo 
objects) with stellar or substellar mass, as well as elementary 
particles (a smooth component). Today, from microlens- 
ing surveys, we have some information about the granular 
component (MACHOs). The absence of very short duration 
events implies that the dark halo cannot be dominated by 
planetary objects. A joint analysis by EROS and MACHO 
collaborations indicated that MACHOs in the mass range 
10" 7 M Q < M < 10" 3 M Q make up less than 25% of the 



dark halo (Alcock et al. 1998). From a likelihood analysis, 
the MACHO collaboration concluded that a population of 
objects of mass ~ 0.5 Mq is consistent with their first two 
year of data. These MACHOs with stellar mass would have 
an important contribution to the total mass (Alcock et al. 
1997; Gould 1997; Sutherland 1999; Mao 2000). However, 
very recent results by the MACHO team, based on approx- 
imately six years of observations, point to a relatively small 
mass fraction (Alcock et al. 2000). For a typical size halo, 
the maximum likelihood estimates suggest the existence of 
a Milky Way dark halo consisting of about 20% MACHOs 
with mass of ~ 0.6 Mq. The EROS collaboration also agrees 
with this small contribution to the dark halo by ~ 0.6 Mq 
objects (Lasserre et al. 2000). Lasserre et al. (2000) derived 
strong upper limits on the abundance of MACHOs with dif- 
ferent masses. For example, < 10% of the dark halo resides 
in planetary objects. Moreover, they ruled out a standard 
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spherical halo in which more than 40% of its mass is made 
of dark stars with 1 Mq . Finally, we remark that the Milky 
Way dark halo inferred from the likelihood method (best 
standard fits by Alcock et al. 2000) is consistent with the 
HST (Hubble Space Telescope) detection of a halo white 
dwarf population (Ibata et al. 1999). A population of cool 
white dwarfs contributing 1/5 of the dark matter in the 
Milky Way could explain all new observational results, but 
the hypothesis presents some difficulties (e.g., Mao 2000; 
Alcock et al. 2000). 

The information on the nature of galaxy dark haloes is 
still largely based on a local spiral galaxy (the Galaxy), and 
so, the study of other galaxies seems an interesting goal. 

The Einstein Cross (QSO 2237+0305) is a z = 1.69 
quasar lensed by a face-on barred Sb galaxy at z = 0.0394. 
The time delay between the four quasar images is expected 
to be less than a day (Rix et al. 1992; Wambsganss & 
Paczynski 1994), and so, one can directly separate intrin- 
sic variability from microlensing signal. For this lens sys- 
tem, light rays of the 4 images pass through the bulge of 
the foreground galaxy and there is a robust evidence that 
microlensing events occur (e.g., Irwin et al. 1989; Wozniak 
et al. 2000). The observed events may be interpreted as a 
phenomenon caused by the granularity of the matter associ- 
ated with the nearby spiral. Very recently, for providing an 
interpretation of the OGLE Q2237+0305 microlensing light 
curve, Wyithe, Turner & Webster (2000) used the contour- 
ing technique of Lewis et al. (1993) and Witt (1993). 

B1600+434 is another interesting gravitational mirage 
lensed by an edge-on disk galaxy. Koopmans & de Bruyn 
(2000) measured the radio time delay between the two im- 
ages of the system and derived a radio difference light curve 
which is in disagreement with zero. They investigated both 
scintillation and microlensing as possible causes of the non- 
intrinsic radio variability. If microlensing is the origin of the 
"anomalous" difference light curve, then it could indicate 
the presence of a lens galaxy dark halo filled with MACHOs 
of mass > 0.5M Q . 

1.2 Microlensing in the first gravitational lens 
system (Q0957+561) 

A third well-known microlensed quasar is the z = 1.41 dou- 
ble system Q0957+561A,B. The main lens galaxy is an ellip- 
tical galaxy (cD) at z — 0.36. While the light associated with 
the image B crosses an internal region of the lens galaxy, the 
light path associated with the component A is « 5 arcsec 
away from the centre of the galaxy. The cD galaxy is close 
to the centre of a galaxy cluster, and consequently, the nor- 
malized surface mass densities ha and kb are the projected 
mass densities of the lensing galaxy plus cluster along the 
lines of sight, normalized by the critical surface mass density. 
Pelt et al. (1998) used the recent values ka = 0.22 and Kb 
— 1.24, which originate from an extended galaxy halo con- 
sisting of the elliptical galaxy halo and an additional matter 
related to the cluster. It is possible that a considerable part 
of the extended halo mass does consist of a dark component, 
although an estimate of the stellar contribution (luminous 
stars) to ka and Kb is not so easy as in the Milky Way. For 
the image B, if the fraction of mass in granular form kbg is 
dominated by normal stars and dark stars similar to the ob- 
jects that have been discovered in the Galaxy (Alcock et al. 



2000), and simultaneously, the main part of the halo mass is 
due to a smooth component {kbg << Kb, kbg << 1) and 
the source quasar is small, then we must expect some long- 
timescale microlensing event caused by one star (luminous or 
dark) crossing the path of this image. In that case of small 
source/one star approximation, the timescale of an event 
will be £ D (years) « 17 ^M(M & ) [600/% (fcms" 1 )], where v t 
is the transverse velocity, and the magnification of the B 
component has a typical duration of several years for a 0.5- 
1 Mq star and any reasonable choice of vt- When kbg is 
high (k B g ~ 1) and/or the source is large, several stars at a 
time must be considered and the model by Chang & Refs- 
dal (1984) is not suitable. The small source/one star model 
by Chang & Refsdal (1984) was generalized in the case of a 
small source and a large optical depth (Paczynski 1986) and 
the case of an extended source and an arbitrary optical depth 
(Kayser, Refsdal & Stabell 1986; Schneider & Weiss 1987; 
Wambsganss 1990). Therefore, the formalisms by Chang & 
Refsdal (1984), Paczynski (1986) and Wambsganss (1990) 
as well as new analytical approximations seem useful tools 
for a detailed analysis of the optical microlensing history of 
QSO 0957+561. A long-timescale microlensing signal was 
unambiguously observed from 1981 to 1999; see Pelt et al. 

(1998) , Press & Rybicki (1998), Serra-Ricart et al. (1999, 
subsequently Paper I). In this paper, we concentrate on the 
possible rapid microlensing signal. In a forthcoming paper, 
we will carry out a comparison between the annual differ- 
ences (averages from January to December) B — A and the 
predictions from different models and physical parameters. 

In the past, using a record of brightness including pho- 
tometric data (in the R band) up to 1995 and a time de- 
lay of 404 days, Schild (1996, hereafter S96) analyzed the 
possible existence of short-timescale microlensing (rapid ex- 
ternal variability on a timescale of months) and very rapid 
microlensing events (with duration of < 3 weeks) in the 
double QSO 0957+561A,B. He found numerous events with 
quarter-year and very short timescales (a few days). S96 also 
claimed that the slower component (events with a width 
of 90 days) can be interpreted as the imprint of an im- 
portant population of microlenses with planetary mass of 
~ 1O~ J M0. Assuming an improved delay value of 417 days, 
Goicoechea et al. (1998, subsequently G98) showed a differ- 
ence light curve corresponding to the 1995/1996 seasons in 
Schild's dataset. G98 obtained fluctuations which could be 
associated with microlensing events, in fact, our results are 
in agreement with the existence of strong microlensing: the 
fluctuations in the difference light curve are clearly inconsis- 
tent with zero and similar to the fluctuations in the quasar 
signal. New work by Schild and collaborators pointed in the 
same direction: adopting a time delay of 416.3 days, Pelt et 
al. (1998) found that Schild's photometry shows evidence in 
favour of the presence of short-timescale microlensing; Schild 

(1999) made a wavelet exploration of the QSO 0957+561 
brightness record, and reported that the rapid brightness 
fluctuations observed in the A and B quasar images (whose 
origin may be some kind of microlensing) are not dominated 
by observational noise; and Colley & Schild (1999), from a 
new reduction of "old" photometric data (subtracting out 
the lens galaxy's light according to the HST luminosity pro- 
file and removing cross talk light from the A and B images 
apertures), derived a structure function for variations in the 
R-band from lags of hours to years, a time delay of 417.4 
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days and a microlensing candidate on a timescale of a day, 
which could imply planetary MACHOs in the lens galaxy 
halo. So, from the photometry taken at Whipple Observa- 
tory 1.2 m telescope by Schild, one obtains two important 
conclusions. First, there is evidence in favor of the exis- 
tence of true short-timescale microlensing signal. Second, 
this rapid signal seems to support the presence of MACHOs 
(in the halo of the cD galaxy) having a very small mass. 
However we note that Gould & Miralda-Escude (1997) have 
introduced an alternative explanation to the possible rapid 
microlensing in the double QSO 0957+561A,B, which is re- 
lated to hot spots or other moving structures in the accre- 
tion disk in the quasar, and so, planetary objects are not 
involved. 

QSO 0957+561A,B was photometrically monitored at 
Apache Point Observatory (Kundic et al. 1995, 1997) in the 
g and r bands, during the 1995 and 1996 seasons. Schmidt & 
Wambsganss (1998, hereafter SW98) analyzed this photom- 
etry and searched for a microlensing signature. Considering 
the photometric data in the g band and a delay of 417 days, 
SW98 produced a difference light curve covering « 160 days 
and concluded that it is consistent with zero. There is no 
variation in the difference light curve with an amplitude in 
excess of ± 0.05 mag and the total magnitude variation of a 
hypothetical microlensing signal is assumed to be less than 
0.05 mag (see the dashed lines in Fig. 1 of SW98). They em- 
ployed this last upper limit to obtain interesting information 
on the parameter pair MACHO-mass/quasar-size. The lack 
of observed fluctuations rules out a population of MACHOs 
with M < 10~ 3 M Q for a quasar size of 10 cm (25%- 100% 
of the matter in compact dark objects). However, other pos- 
sible scenarios (e.g., a small source and a halo consisting of 
MACHOs with M > 10~ 2 M , a source size of 10 15 cm and 
a halo with compact dark objects of mass < 10~ 3 M , etc.) 
cannot be ruled out from the bound on the microlensing 
variability in the 160 days difference light curve. In short, 
SW98 have not found reliable evidence for the presence of 
rapid microlensing events. 

The gravitational lens system Q0957+561 was also 
monitored by our group with the IAC-80 telescope at 
Teide Observatory, from the beginning of 1996 February to 
1998 July (see Paper I). We re-reduced our first 3 seasons 
(1996-1998) of QSO 0957+561 observations in the R band, 
made the difference light curves for 1996/1997 seasons and 
1997/1998 seasons and studied the origin of the deviation be- 
tween the light curves of the two images. All the results are 
presented in this article. The plan of the paper is as follows: 
in Sect. 2 we present the difference light curves and report 
on new constraints on microlensing variability. In Sect. 3 we 
suggest different models that explain the difference signal. 
In Sect. 4 we discuss the sensitivity of the telescope to dif- 
ferent microlensing "peaks". At the end of the paper (Sect. 
5), we summarize our results. 



2 FIRST THREE SEASONS OF QSO 0957+561 
OBSERVATIONS IN THE R BAND: 
DIFFERENCE LIGHT CURVES 

We have been monitoring Q0957+561 over the past 4 years 
(from 1996 February) with the 82 cm IAC-80 telescope 
(at Teide Observatory, Instituto de Astrofisica de Canarias, 



Spain) and have obtained a large R band dataset. The con- 
tribution to the solution of the old controversy regarding the 
value of the time delay (« 400-440 days or > 500 days ?) 
was the first success of the monitoring program (Oscoz et al. 
1996; see also Kundic et al. 1995, 1997; Oscoz et al. 1997). 

In order to give refined measurements of both time 
delay and optical microlensing, we have introduced some 
modification with respect to the original aperture photom- 
etry (see Oscoz et al. 1996). Reduction of the images A 
and B is complicated by the presence of cross contami- 
nation and contamination from light of the main lensing 
galaxy. The two kinds of contamination depend on the see- 
ing, and it is not clear what is the optimal way of ob- 
taining the best photometric accuracy. At present, we re- 
duce each available night by fitting a profile to the im- 
ages, which is consistent with the point spread function of 
comparison stars. This new method of reduction and the 
photometry from 1996 to 1998 (the first 3 seasons) are de- 
tailed in Paper I. A table including all data is available at 
http : //www . iac . es/pro j ect/quasar/lens7 . html. 

In the QSO 0957+561 quasar, a time delay of w 420 
days is strongly supported (e.g., G98). Using our first 3 sea- 
sons of data, the time delay estimates (in Paper I) are of 425 
± 4 days (from the <5 2 -test, which is based on discrete corre- 
lation functions) and 426 ± 12 days (from dispersion spec- 
tra). A comparison between the discrete cross-correlation 
function and the discrete autocorrelation function, indicates 
that a time delay of < 417 days is in disagreement with the 
photometry (see Fig. 16 of Paper I) , while a delay of about 
425 days is favoured. Thus, we adopted a time delay of 425 
days. 

We concentrate now on the difference light curves. In 
order to estimate the difference light curve (DLC) for the 
1996/1997 and the 1997/1998 seasons, we used 30 observa- 
tions of image A corresponding to the 1996 season (A96), 
28 observations of image B corresponding to the 1997 season 
(B97), 44 photometric data of image A in the 1997 season 
dataset (A97) and 84 photometric data of image B in the 
1998 season dataset (-B98). There are about 100 days of over- 
lap between the time-shifted (in the time delay) light curve 
A96 and the light curve B97, and about 230 days of overlap 
between the time delay-corrected light curves A97 and B98. 
A main problem of the IAC-80 telescope (using the avail- 
able observational time of 20-30 min/night) is related to the 
photometric errors. The mean errors in the initially selected 
datasets are approximately 19 mmag (^496), 24 mmag (B97), 
28 mmag (A97) and 24 mmag (B98). For short-timescale 
microlensing studies, these errors are large and one must re- 
reduce the data (grouping them for obtaining lower errors) . 
Because of the possible rapid microlensing variability on one 
month timescale, the timescale of the groups should not be 
too large (< 10 days); it should not be too small for having 
a sufficient number of data, and so, relatively small errors. 
The re-reduced photometry consists of 12, 11, 22 and 36 "ob- 
servations" in four new (and final) datasets A96, B97, A97 
and B98, respectively. These datasets are available by send- 
ing a request to rmerino@astro.physik.uni-potsdam.de. 
For groups in A96, the timescales are less than 3 days and 
the mean error is of ~ 12 mmag, for groups in B97, the 
timescales are < 8 days and the mean uncertainty is of ~ 16 
mmag, for groups in A97, the timescales are also < 8 days 
and the mean error is lowered to £s 20 mmag, and for groups 
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Figure 1. Difference light curve for 1996/1997 seasons (in the R 
band) . We used bins with semi-size of 9 days and adopted a time 
delay of 425 days. The times associated with the circles are the 
dates in the time-delay shifted light curve A96 (see main text). 

in B98, the maximum timescale and the mean uncertainty 
are 6 days and 16 mmag, respectively. Therefore, making 
groups with a maximum timescale of w 1 week (the mean 
timescale is of ~ 2 days), the mean errors are lowered in 7-8 
mmag. 

As we have seen in the previous paragraph, the bright- 
ness record for the A image (A96 or A97) is measured only 
at a set of discrete times U (i = 1,...,N) and the light curve 
for the B image (B97 or B98) is also determined at dis- 
crete times tj (j — 1,...,M). Since the observational light 
curves are irregularly sampled signals, to obtain the DLC 
(A96/B97 or A97/B98), we can use different methodolo- 
gies, for example, the interpolation suggested by SW98 or 
the binning that appears in G98. Here, we are interested in 
the DLC binned in intervals with size 2a around the dates 
in the light curve A TS (time delay-shifted light curve A). In 
other words, each photometric measurement Aj s at the date 
ti + Atba, where Atba is the time delay, will be compared 
to the observational data Bf s = Bj + <A>-<B><i,tti 
+ Atba - a < tj < U + Atba + ce (B MS is the magnitude- 
shifted light curve B). The values Bj 1 " 5 within each bin are 
averaged to give < Bf IS >t (i = 1,...,N), and one obtains 
the difference light curve (DLC) 

Si =< Bf s > t -Aj s , (1) 

being i = 1,...,N. The observational process A TS (t) can be 
expanded as an intrinsic signal s(t) plus a noise process 
riA(t) related to the procedure to obtain the measurements, 
and a microlensing signal 7?ia(£)- In a similar way, B MS (t) 
= s(t) + iis(t) + ms(t). So, the deviation Si must be inter- 
preted as a combination of several factors, i.e., 

Si = [< sj >i — Si] + [< n B j >i — nAi] + [< m B j >i — mAi]-(2) 

If s(t) is a smooth function, then s< = s(ti) and Sj = s(tj), 
while when s(t) is a stochastic process, Si represents a re- 
alization of the random variable s(U) and Sj denotes a re- 
alization of the random variable s(tj). With respect to the 
observational noise, ua% is a realization of the random vari- 
able nA{U) [similarly, nsj is one of the possible values of 
ns(ij)]. Also, in Eq. (2), rriAi = mA(U) and msj = msfe). 
From Eq. (2) it is inferred that the difference signal will be 



Figure 2. Difference light curve for 1997/1998 seasons (in the R 
band) . We used bins with semi-size of 8 days and adopted a time 
delay of 425 days. The deviations [5;; see Eq. (1)] are evaluated at 
discrete dates corresponding to the time-delay shifted light curve 
A97. 

never zero, even in absence of microlensing. There is a back- 
ground dominated by observational noise, which is present 
in any realistic situation. In the case of very weak or null mi- 
crolensing, we expect a trend of the DLC rather consistent 
with zero [taking into account the standard errors ei,...,£jv 
in the deviations estimated from Eq. (1)]. However, in the 
case of strong microlensing, several absolute deviations |<5i| 
should be noticeably larger than the associated uncertainties 
£i- 

For the 1996/1997 seasons (from the final datasets A96 
and B97), using a time delay of 425 days and bins with semi- 
size of a = 9 days, we derived the DLC that appears in Fig. 
1. Two thresholds are also illustrated: ± 0.05 mag (discon- 
tinuous lines). In Fig. 1, there is a "peak" around day 1615: 
two contiguous points significantly deviated from the zero 
line, that verify \Si\ > e;. If the whole DLC is modelled as a 
single Gaussian event and the data are fitted to the model, 
we obtain that the amplitude and the full-width at one- 
tenth maximum (FWTM) of the Gaussian law must be w - 
33 mmag and w 14 days, respectively (best-fit characterized 
by X 2 /N ~ 1). Apart from this very short duration event, 
which is probably caused by observational noise (see next 
section) , there is no evidence in favor of the existence of an 
event on longer timescales. We note that "event" is used in a 
general sense, and it may be due to true microlensing, obser- 
vational noise, a combination of both or other mechanisms. 
In particular, none Schild-event (events having a width of 
three months and an amplitude of ± 50 mmag; see S96) is 
found. Although the difference signal is only tested during a 
100 days period, to find a Schild-event belonging to a dense 
network of similar fluctuations (positive and negative), the 
"sampling" would be sufficient. In any case, from our second 
DLC (see here below), we must be able to confirm/reject the 
existence of a network of events with quarter- year timescale 
and an amplitude of ± 50 mmag. Finally, there are derived 
bounds on the amplitude of the microlensing fluctuations of 
± 0.05 mag, which are similar to the bounds for 1995/1996 
seasons (see Introduction and SW98). 

For the 1997/1998 seasons (from the final datasets A97 
and B98), we also made the corresponding DLC. In Fig. 2, 
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the DLC and two relevant thresholds are depicted. The dif- 
ference signal is in apparent agreement with zero, i.e., Fig. 
2 shows a noisy relationship B MS = A TS . We observe no 
Schild-events, and therefore, the total difference signal (ml 
year of overlap between the time-shifted light curve for the 
A component and the magnitude-shifted light curve for the 
B component) is in clear disagreement with the claim that 
90 days and ± 50 mmag fluctuations occur almost contin- 
uously. One can also infer constraints on the microlensing 
variability. In good agreement with the DLCs for 1995-1997 
seasons, a hypothetical microlensing signal cannot reach val- 
ues out of the very conservative interval [- 0.05 mag, + 0.05 
mag] . We finally remark that the methodology introduced by 
SW98 (the technique of interpolation) leads to DLCs similar 
to the DLCs discussed here (Figs. 1-2) . 



3 INTERPRETATION OF THE DIFFERENCE 
SIGNAL 

The DLCs presented in Sect. 2 are in apparent agreement 
with the absence of microlensing signal. However, to settle 
some doubts on the ability of the observational noise in order 
to generate the observational features (e.g., the very rapid 
event in Fig. 1) and the measured variabilities (rms aver- 
ages), a more detailed analysis is needed. In this section, we 
are going to test three simple pictures without microlensing. 
In brief, the ability of some models for generating combined 
photometries and difference signals similar to the observa- 
tional ones is discussed in detail. 

The observational combined photometry consists of 
both light curves A TS and B MS . Thus, assuming that m(t) 
= 0, the combined light curve (CLC) must be related to a 
process C(t) = s(t) + n(t). The intrinsic signal s(t) is cho- 
sen to be either a smooth function (polynomial; picture I) 
or a polynomial plus a stationary noise process (picture II) 
or a correlated stationary process (picture III). In the first 
case (picture I), we work with s(t) = X/p-o a P^ P ( wnen the 
CLC is reasonably smooth, this intrinsic signal is a suitable 
choice). The polynomial law leads to Ck = J^"_ a p t p k + 
rife at a date tk, where Ck (k = 1,...,N+M) are the com- 
bined photometric data. Considering that the process n(t) 
is Gaussian with < n(t) > = and a„{t) = < n 2 (t) >, 
and identifying the measurement errors a\ with the noise 
process cr 2 (tk), the probability distribution of nk at a given 
time tk is Pk{n k ,t k ) = (\j \THok) exp(-n|/2<Tfe). Here, the 
angle brackets denote statistical expectation values. As the 
random variables n(tk), k = 1,...,N+M, are independent (the 
noise is uncorrelated with itself), the joint probability dis- 
tribution of the noise vector n = (m,...,riN+M) is given by 

N+M 

P(n) = J] P k 

k=l 

N+M n 

= (2^^ H(l/a k )eM-[Ck-J2 a A] 2 /^t}- (3) 

k=l p=0 

Maximizing the likelihood function L =lnP with respect 
to the parameters a p , or equivalently, minimizing \ 2 = 
Ylk-i^jCk — Xr=o a p^fe] 2 /°fc> we nn d a possible reconstruc- 
tion of the intrinsic signal (and thus, a model). If this pro- 
cedure does not work (e.g., \ 2 /dof is relatively large, with 



dof = N+M-(n+l) being the number of degrees of freedom), 
we perform a fit including a stationary intrinsic noise as an 
additional ingredient (picture II). This new ingredient can 
account for noisy CLCs. The intrinsic noise n(t) is taken 
to be Gaussian with < n{t) > — and a^(t) = a 2 nt , and 
moreover, n(t) is uncorrelated with both n(t) and with it- 
self. Now, C(t) = s(t) + £(t), where s(t) = YZ=o a P tP and 
£(t) = n(t) + rj(t), and we focus on the global noise pro- 
cess £(t). As the processes n(t) and r/(t) are Gaussian and 
mutually independent, their sum is again Gaussian, and the 
average and variance of £ (t) are the sums of the averages and 
variances of both individual noise processes. The probability 
distribution of £k at an epoch tk can be written as Pk(^k, tk) 
= [lA/2^(ai+aL) 1/2 ] exp[-^/2(^+cTL)], and the joint 
probability distribution of the noise vector £ = (£i,...,£jv+m) 
should be P(£) — Ylk^ Pk(£k,tk)- Finally, instead of the 
standard procedure (to maximize the likelihood function), 
we equivalently minimize the function 

N + M n 

X = ^2 {ln(°fc + Vint) + \ C k - ^ a P tF k\ 2 / ( a k + a int)}-{^) 
k=l p=0 

Through this method, the intrinsic signal is partially recon- 
structed. We find the coefficients of the polynomial and the 
variance of the intrinsic noise, but after the fit, the realiza- 
tions r]k (k — 1,...,N+M) remain unknown. However, the de- 
rived model permits us to make simulated CLCs and DLCs, 
since only the knowledge of the smooth intrinsic law and the 
statistical properties of the noise processes are required for 
this purpose. 

A very different procedure was suggested by Press, Ry- 
bicki & Hewitt (1992 a,b, hereafter PRH92). They assumed 
the intrinsic signal as a correlated stationary process. For 
this case III, it is possible a reconstruction of the realiza- 
tions of s(t), provided that the correlation properties are 
known. PRH92 considered that the observational noise n(t) 
is uncorrelated with s(t) (and with itself), and therefore, 
only the autocorrelation function K s (t) = < s(t)s(t + r) > 
is needed, being s(t) = s(t)— < s >. The autocorrelation 
function of the intrinsic signal is not known a priori and 
must be estimated through the CLC. We can relate the au- 
tocorrelation properties to the first-order structure function 
D«(r)by 

fl( 1 »(T) = (l/2 1 -)^( Sm - S( ) 2 

l,m 

« \ < W + r)- s(t)] 2 >= K s (0) - K s (r), (5) 

where the sum only includes the (l,m) pairs verifying that 
t m — U w r (the number of such pairs is v) . From the CLC, 
one infers (e.g., Haarsma et al. 1997) 

D«(r) « (1/21/) ^[(C™ - - "I - ° 2 ™\> ( 6 ) 

l ,m 

which is an evaluation of the difference K a (0) - K e (r). As 
usual we assume a power-law form for the first-order struc- 
ture function, and perform a fit to the power law. Finally, 
the variance of the intrinsic process K s (0) is assumed to be 
the difference between the variance of the CLC and a cor- 
rection due to the observational noise. The whole technique 
is described in PRH92 and other more recent papers (e.g., 
Haarsma et al. 1997). 
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Figure 3. The combined photometry of QSO 0957+561A,B for 
the 1996/1997 seasons in the R band (at Teide Observatory). 
The open circles trace the time-shifted (+ 425 days) light curve 
A96 and the filled squares trace the magnitude-shifted (+ 0.0658 
mag) light curve B97. The lines are related to two reconstruc- 
tions of the intrinsic signal: considering an intrinsic signal of the 
kind polynomial plus stationary noise (top panel) and the optimal 
reconstruction following the PRH92 method (bottom panel). 

3.1 The 1996/1997 seasons 

For the 1996/1997 seasons, we first have done the corre- 
sponding CLC. In a second step, using the picture I (see 
above), we attempted to fit the combined photometry. A 
quadratic law (n — 2) gives \ 2 /dof = 1.65 (best fit), whereas 
X 2 /dof(n = 1) = 2.52, X 2 /dof(n = 3) = 1.74 and X 2 /dof(n 
= 4) = 1.83. Thus the modelling of the CLC has proven 
to be some difficult. Fortunately, the inclusion of an intrin- 
sic noise (picture II) with moderate variance does not fail 
to generate an acceptable fit. When the intrinsic signal is 
the previous best quadratic fit to which an intrinsic noise 
with <j int = 9 mmag is added, we obtain X 2 /dof = 1.15 
(x 2 /N+M = 0.95). The quality of the fit has changed signif- 
icantly with the addition of the new noise, whose variance 
(Tint = 9 mmag) is less than the mean variance of the ob- 
servational noise (= 12-16 mmag). In Fig. 3 (top panel) the 
CLC and the reconstruction are presented. The open circles 
represent the time-shifted light curve A96, while the filled 
squares are the magnitude-shifted light curve B97. The best 
polynomial (n = 2) is traced by means of a solid line, and 
the two lines with points are drawn at ± 9 mmag (the best 
value of o in t) from the polynomial. Apart from the CLC, 
we checked the observational structure functions [see 
Eq. (6)] and D?' as well as the predictions (with respect 
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Figure 4. The first-order and second-order structure functions 
(1996/1997 seasons in the R band). The open circles are the values 
inferred from the observational data and the filled triangles are 
the predictions from the reconstruction of the kind polynomial + 
stationary noise. The observational first-order structure function 
was fitted to a power-law Et e (solid line in the top panel). As- 
suming this fit as an estimation of the autocorrelation properties 
of a hypothetical correlated stationary process (K s (0) - K s (t)), 
the predicted second-order structure function is illustrated by a 
solid line in the bottom panel. 

to the structure functions) from our first succesful recon- 
struction. The observational second-order structure function 
is computed in the following way (see Simonetti, Cordes & 
Heeschen 1985; we take a normalization factor equal to 1/6): 

D?\t) « (1/6/i) [(C*n-2C m + C i ) 2 -of -4a 2 m -a 2 n ],(7) 

where fx is the number of (l,m,n) valid triads so that t m — tj ~ 
r and t„ — U « 2r. Also, the predicted structure functions 
are 

dP(t) » (l/2i/) 5}«(t m ) - S(tj)] a + °lnU 

l.m 

D {2) {t) » (1/6/*) _ 2S (*"*) + S (*')] a + ( 8 ) 

being s(t) the fitted quadratic law. Fig. 4 shows the good 
agreement between the observational values (open circles) 
and the predicted trends (filled triangles). This result con- 
firms that the reconstruction is reliable. The meaning of the 
two straight lines in Fig. 4 will be explained here below. 

Our interest in this paper is less directly in the de- 
tails of a given reconstruction of the underlying intrinsic 
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Figure 5. Global properties of the measured photometry for 
1996/1997 seasons (filled star) and 100 simulated photometries 
(open circles). The numerical simulations arise from Ml, which is 
a model with three ingredients: polynomial law + intrinsic noise 
+ observational noise. 



Figure 7. Gaussian events (they are classified according to their 
amplitude and FWTM) found in the first 33 simulations via Ml. 
The number of features as well as the amplitudes and time-scales 
are relatively surprising. 
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Figure 6. The true DLC for 1996/1997 seasons (left-hand top 
panel) together with 3 simulated DLCs (via Ml). The solid lines 
are fits to Gaussian events. A curious result observed in the sim- 
ulated DLCs is the existence of events, which could be naively 
interpreted as microlensing fluctuations. 



signal than it is in analyzing simulated photometries con- 
sistent with the reconstruction and with the same sampling 
(dates) and errors as the measured data. The first model 
(Ml) comprises the best quadratic fit in the absence of in- 
trinsic noise (a smooth component) and a Gaussian noise 
process characterized by a known variance at discrete times 
tk- <y\ + (Ti nt . From Ml we derived 100 simulated CLCs 
and the corresponding DLCs. We remark that, in each sim- 
ulation (CLC), N simulated data points represent a syn- 
thetic light curve A TS , while the other M data are simu- 
lated measurements of the magnitude-shifted light curve B. 
Fig. 5 shows the relationship between the values of x 2 /N+M 
(X 2 = J2k=i M \- Ck ~ ^(ife)] 2 /If* + &int]) and the rms aver- 
ages of the DLCs (rms = [i J2?=i S i] 1/2 )- The 100 °P en 
circles are associated with the simulated photometries and 
the filled star is related to the measured photometry. The 
true (measured) photometry appears as a typical result of 



the model. One sees in the figure a broad range for CLC- 
X 2 /N+M (0.2-2.2) and DLC-rms (8-36 mmag), and the true 
values of CLC-x 2 /N+M = 0.95 and DLC-rms = 22 mmag 
are well placed close to the centre of the open circle distri- 
bution. Thus, the measured combined photometry seems a 
natural consequence of Ml, which is a model without very 
rapid and rapid microlensing. However, due to the event 
found in Fig. 1 (around day 1615) and other local features 
less prominent than the event, we would be to doubt this 
conclusion and to study details in the synthetic DLCs. In 
Fig. 5, to provide some guidance, the open circles corre- 
sponding to simulated datasets with CLC-x 2 /N+M similar 
to the measured value have been enclosed in a rectangular 
box. Also, we have drawn an elliptical surface centred on the 
filled star, which includes (totally or partially) three open 
circles associated with the synthetic photometries analogous 
(global properties of both the CLC and the DLC) to the true 
brightness record. As we must put into perspective the very 
rapid event and other local properties discovered in the true 
DLC for 1996/1997 seasons, this DLC and its features were 
compared with the three DLCs that arise from the simula- 
tions. In Fig. 6 we present the comparison. All events (each 
event includes a set of two or more consecutive deviations 
which have equal sign and are not consistent with zero) has 
been fitted to a Gaussian law and marked in the figure. The 
measured DLC (left-hand top panel and Fig. 1) is not differ- 
ent to the other three. In fact, in our 100 th simulation (slOO; 
right-hand bottom panel), two events appear. The positive 
event is more prominent than the negative event, and this 
last one is similar to the measured one. With respect to 
the regions without events, the true variability cannot be 
distinguished from the simulated ones. To throw more light 
upon the problem, we searched for Gaussian events in 1/3 
of all simulations (sl-s33), as a sample of the whole set of 
simulations because the computation turned out to be very 
time-consuming. The results are plotted in Fig. 7: ampli- 
tude of each event (mmag) vs. FWTM (days). There are a 
lot of events with amplitude in the interval [- 50 mmag, + 
90 mmag] and duration < 70 days. In particular, the proba- 
bility of observing a negative event is of 15% and the prob- 
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Figure 8. Global properties of the true DLC for 1996/1997 sea- 
sons (filled star) and 100 simulated DLCs (open circles). The 
numerical simulations were made through a model including the 
optimal reconstruction of a correlated stationary process and a 
Gaussian observational noise process whose variance at the dates 
of the real data is known. 
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Figure 9. Four simulated DLCs (via M2). For comparison with 
the true event in Fig. 1 (see also Fig. 6), the Gaussian events have 
been clearly marked on the panels. 



ability of observing one or more events is of about 50%. So, 
it must be concluded that the noisy (around zero) difference 
light curve based on observations is totally consistent with 
Ml, and the deviations from the zero line can be caused by 
the combined effect of the processes n(t) (main contribution) 
and rj(t). 

We also propose a reconstruction of the underlying in- 
trinsic signal as realizations of a correlated stationary pro- 
cess (picture III). The observational first-order structure 
function can be fitted to a power-law Er e (see Fig. 4, solid 
line in the top panel) . If one considers this fit as an estimate 
of the difference K s (0) - K s (t) for a correlated stationary 
process, then it is straightforward to obtain the predicted 
second-order structure function (see Fig. 4, solid line in the 
bottom panel: the prediction is irrelevant to reconstruct the 
intrinsic signal, but it is necessary for testing the consis- 
tency of the starting point K s (0) - K s (t) — Et e ) and to 
apply the reconstruction formalism by PRH92. Therefore, 



Figure 10. Gaussian events found in the numerical simulations 
(via M2) with 20 < rms (mmag) < 24 (open circles). Events 
very similar to the real event (filled star) are produced in two 
simulations. 



we are able to find the realizations of the intrinsic process 
at the observational times t k (k — 1,...,N+M) as well as in 
the gaps between the observations. The PRH92 technique 
leads to an acceptable fit with x 2 /dof — 1.18 (dof = N+M- 
1), and our second successful reconstruction is showed in 
Fig. 3 (bottom panel). The knowledge of both the optimal 
reconstruction and the properties of the Gaussian observa- 
tional noise process at discrete times tk (k — 1,...,N+M), 
permits us to make 100 new simulations. In Fig. 8 details of 
the rms averages of the DLCs are provided (open circles). 
The observational DLC has a rms average (filled star) simi- 
lar to the rms average of about 1/5 (20%) of the simulated 
DLCs. Furthermore, four simulated DLCs with rms in the 
interval [20 mmag, 24 mmag] (in Fig. 8, this range of vari- 
ability is labeled with two horizontal lines) appear in Fig. 
9. From the new model (M2), DLCs with no events (as in 
the analysis presented above, the Gaussian events are re- 
lated to "peaks", or in other words, we only made events 
around consecutive multiple deviations with equal sign and 
well separate from zero) and DLCs that incorporate more or 
less prominent features are derived. We note that one DLC 
(right-hand bottom panel) has an event almost identical to 
the true one in Fig. 1 and Fig. 6. Fig. 10 shows the proper- 
ties of all Gaussian events in the simulated DLCs with rms 
in the vicinity of the observational rms (open circles). The 
measured event is also depicted (filled star), and we can see 
two simulated events analogous to it. We finally conclude 
that the observational DLC is in clear agreement with M2, 
and so, microlensing would be not advocated. In this frame- 
work (M2), the observational noise process is a sufficient 
mechanism for originating the measured deviations. 



3.2 The 1997/1998 seasons 

The combined photometry for 1997/1998 seasons and the re- 
construction based on a polynomial fit are showed in Fig. 11 
(top panel). The open circles represent the time-shifted light 
curve ^497 and the filled squares are the magnitude-shifted 
brightness record B98. There is no need for the presence of 
an intrinsic noise, and a simple quadratic law works well, 
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Figure 11. The combined photometry for 1997/1998 seasons in 
the R band (at Tcidc Observatory). The open circles trace the 
time-shifted (+ 425 days) light curve A97 and the filled squares 
trace the magnitude-shifted (+ 0.0603 mag) light curve B98. The 
solid lines represent two reconstructions of the intrinsic signal: 
the best quadratic fit (top panel) and the optimal reconstruction 
following the PRH92 method (bottom panel). 



leading to x 2 /dof = 0.85 (best fit). In Fig. 11 (top panel), 
the solid line traces the reconstruction of the intrinsic signal. 
Besides the comparison between the measured CLC and the 
fitted polynomial, we tested the predicted structure func- 
tions. In Fig. 12 we present the observational and 
[open circles; see Eqs. (6-7)] and the predictions from the 
best quadratic fit (filled triangles; see Eqs. (8) with a int = 
0) . The laws traced by the dashed and solid lines in this fig- 
ure will be discussed below. It is evident that the behaviours 
deduced from observations and the predicted trends agree 
very well, and this result indicates that the reconstruction 
is robust. 

The first model for 1997/1998 seasons (M3) consists of 
the best quadratic fit together with a Gaussian observational 
noise process (whose variance is a\ at discrete times t k , k 
= 1,...,N+M, being at the measurement errors at the dates 
of observation tk)- Using M3 we performed 100 simulated 
CLCs (and consequently, 100 simulated DLCs). The global 
properties of the simulated photometries (open circles) and 
the true dataset (filled star) are depicted in Fig. 13. If we 
concentrate on the simulations with x /N+M similar to the 
measured value (rectangular box), the true DLC has a rms 
relatively small (of about 15 mmag), but consistent with 
the rms distribution associated with the simulated DLCs. 
We remark that 3 simulations (open circles in the ellipti- 
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Figure 12. The first-order and second-order structure functions 
(1997/1998 seasons in the R band). The open circles are the val- 
ues inferred from the observational data and the filled triangles 
are the predictions from the reconstruction of the kind polyno- 
mial. The observational first-order structure function was fitted 
to different laws, and two "reasonable" fits are drawn in the top 
panel (dashed and solid lines). If the fits arc interpreted as the 
difference K s (0) - K s (t) for a correlated stationary process, the 
corresponding predicted second-order structure functions are il- 
lustrated by two lines in the bottom panel. 
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Figure 13. Global properties of the true photometry for 
1997/1998 seasons (filled star) and 100 simulated photometries 
(open circles). The numerical simulations are based on a model 
of the kind polynomial plus observational noise. 
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Figure 14. The true DLC for 1997/1998 seasons (left-hand top 
panel) together with 3 simulated DLCs (via M3). An only "peak" 
is marked by a double arrow (see right-hand top panel). 



Figure 15. Global properties of the true DLC for 1997/1998 
seasons (filled star) and 100 simulated DLCs (open circles). The 
numerical simulations were made from M4 (see main text). 



cal surface around the filled star) are analogous to the real 
brightness record, and in Fig. 14, their DLCs can be com- 
pared with the true DLC. The measured difference signal 
(left-hand top panel and Fig. 2) is a quasi-featureless trend 
and similar to the other synthetic DLCs. There are no signifi- 
cant events in these four DLCs with small global variability. 
We conclude that a model with no microlensing (M3) has 
the ability of generating ligth curves like the real data for 
1997/1998 seasons. Henceforth, we are going to treat the 
"peaks" as top-hat fluctuations, i.e., given a "peak" includ- 
ing deviations Spi,...,6pp at times tpi,...,tpp, the amplitude 
and duration of the associated top-hat profile will be evalu- 
ated as the average of the individual deviations and the dif- 
ference tpp — tpi, respectively. In Fig. 14, a "peak" (defined 
by two contiguous negative deviations, which are inconsis- 
tent with zero) appears in the DLC from the 7 th simulation 
(s7; right-hand top panel). The "peak" is marked by a dou- 
ble arrow that represents the amplitude and duration of the 
associated top-hat profile. 

The inspection of the observational first-order struc- 
ture function (see Fig. 12) suggests that the underlaying 
law could be intricate. To find the autocorrelation proper- 
ties of a possible and plausible correlated stationary process 
causing the main part of the observed signal (picture III), 
this observational structure function was firstly fitted to a 
non-standard law (t) = Et s /[1 + (t/T) x ] 2 . As showed 
in Fig. 12 (dashed line in the top panel), the fit is excellent. 
However, when we attempt to reproduce the observational 
second-order structure function, an inconsistent prediction 
is derived (dashed line in the bottom panel). The predic- 
tion fails at r < 70 days. Other functions led to fits more 
or less successful, and finally we adopted the point of view 
by PRH92. In Fig. 12 (top panel) one sees a power-law be- 
haviour up to r = 140 days. The drop at the largest lags is 
due to the coincidence of values in the starting and ending 
parts of the measured CLC. Therefore, we assume that the 
observational first-order structure function is a reliable es- 
timator of K a (0) - K s (t) at t < 140 days, whereas it is a 
biased estimator at r > 140 days. The power-law fit to the 
data at lags r < 140 days gives the autocorrelation proper- 
ties for the correlated stationary process, shown as a solid 




Figure 16. Four simulated DLCs via M4. No events are found 
(for a comparison, see Fig. 2). 

line in the Fig. 12 (top panel). The predicted second-order 
structure function (Fig. 12, solid line in the bottom panel) 
is consistent with the observational one up to a lag of 70 
days, and it deviates from the observational trend at r > 70 
days. However, since the observational second-order struc- 
ture function at lag r is associated with the autocorrelation 

(2) 

at lag 2t, the observational Dl ' (r > 70 days) will be re- 
lated to the autocorrelation at r > 140 days, which is poorly 
traced from observations. Thus the deviation at largest lags 
is reasonable and the global prediction should be considered 
as a consistent result. 

Once the relationship between the structure function 
and the autocorrelation has been established, we can directly 
obtain both an optimal reconstruction of the realizations of 
the intrinsic signal and a new model (M4). The relatively 
smooth reconstruction is showed in Fig. 11 (bottom panel; 
the x 2 /dof value is of 0.86), and the associated model leads 
to 100 simulations, whose global properties (rms averages of 
the DLCs) are presented in Fig. 15 (open circles). In Fig. 
15, a filled star represents the true rms average, which is 
consistent (although marginally) with the rms distribution 
from simulations. Finally, four simulated DLCs with rms < 
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Figure 17. Probability distributions of the rms averages of the 
synthetic DLCs. The numerical simulations were made from Ml 
(solid line) and M2 (dashed line). 
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Figure 19. Top-hat fluctuations found in the numerical simu- 
lations based on Ml. We show 84 features that appear in 100 
simulated DLCs. 
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Figure 18. Probability distributions of the rms averages of the 
synthetic DLCs. The numerical simulations are based on M3 
(solid line) and M4 (dashed line). 
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Figure 20. Top- hat fluctuations in 100 simulated (via M2) DLCs. 
They were found 55 " peaks" . 



17 mmag (in Fig. 15, the upper limit of 17 mmag is marked 
with one horizontal line) have been selected for a more de- 
tailed inspection. We found noisy behaviours around zero 
and no events in these synthetic DLCs, i.e., the results agree 
with the analysis of the real difference signal for 1997/1998 
seasons. The 4 quasi-featureless simulated DLCs appear in 
Fig. 16. We again showed that microlensing is not necessary. 
The real combined photometry and difference signal can be 
due to a set of realizations of two very different processes: 
a correlated stationary process (intrinsic) and a Gaussian 
noise (observational). 



4 THE ABILITY OF THE IAC-80 TELESCOPE 
TO DETECT MICROLENSING "PEAKS" 

The sensitivity of our telescope to microlensing variabil- 
ity in a given observational DLC is an important issue 
which merits more attention. To explain the observations for 
1996/1997 seasons and 1997/1998 seasons, we proposed (in 
Sect. 3) four models based on pictures including only an in- 



trinsic signal and observational noise. The simulations aris- 
ing from these models (100 simulated difference light curves 
per model) are a useful tool to study the statistical proper- 
ties of the expected difference signal in the absence of mi- 
crolensing, and so, to test the resolution of the IAC-80 tele- 
scope for microlensing variability. In Fig. 17 we present the 
probability distributions of the rms values (DLCs) derived 
from Ml (solid line) and M2 (dashed line). A value of about 
20 mmag has a relatively high probability of 20-40%, while 
a rms exceeding 36 mmag is inconsistent with both models, 
as can be seen in Fig. 17. Fig. 18 also shows the probability 
of observing (in the absence of microlensing) different rms 
values: via M3 (solid line) and via M4 (dashed line). The rms 
averages in the interval 19-27 mmag are highly probable (20- 
40%), but a global variability characterized by either rms < 
12 mmag or rms > 38 mmag can be excluded. As a general 
conclusion, the rms of the difference signal induced by noise 
does not exceed a threshold of 37 mmag. Therefore, the rms 
values of future observational DLCs can be used to discrimi- 
nate between the presence of the expected background (global 
variability with rms < 40 mmag) and the probable existence 
of true microlensing signal (rms > 40 mmag). 
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Figure 21. Top-hat fluctuations from M3. We note the existence 
of noise "peaks" with a duration longer than 40 days. All these 
features are however associated with an unfortunate small gap in 
our photometry. 
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Figure 22. "Peaks" from M4. 

The previous discussion on the global variability is in- 
teresting, but it is not the main goal of this section. Our 
main goal lies in discussing the sensitivity of the telescope 
(taking into account typical sampling, photometric errors, 
re-reduction and making of bins) to several classes of mi- 
crolensing "peaks" (the cores of the microlensing events). 
We have seen, in Sect. 3, a figure that shows the properties 
of the Gaussian events (amplitude and FWTM) found in a 
subset of simulations from Ml (see Fig. 7). Fig. 7 can be 
compared to the distribution of top-hat fluctuations found 
in all DLCs generated with Ml. In Fig. 19 the distribu- 
tion of the top-hat fluctuations (basically the properties of 
the "peaks" associated with them) appears, and a direct 
comparison between Fig. 7 and Fig. 19 indicates the log- 
ical fact that Gaussian fits lead to longer durations than 
top-hat estimates. In the case of Gaussian fits, events with 
a duration (FWTM) of 1-2 months are abundant and only 
features with a timescale > 70 days are ruled out. However, 
the "peaks" (from Ml) with a timescale of about one month 
are scarce. To discuss the power of resolution of the telescope 
for local microlensing variability we chose the top-hat fluc- 
tuations ("peaks") instead of the events. The properties of 



an event (around a "peak") depend on the assumed profile 
(e.g., Gaussian, Lorentzian, etc.) and the global behaviour 
of the DLC, whereas the top-hat shape directly traces the 
"peaks", avoiding to make assumptions on their wings and 
the use of the rest of the corresponding DLCs. In a few 
words, the top-hat structures are more local and free from 
assumptions than the events. 

The "peaks" from M2 (Fig. 20) are not so numerous as 
the top-hat fluctuations inferred from the first model (Ml). 
Moreover, the new cloud of points (open circles) is more 
concentrated towards shorter durations. In fact, all "peaks" 
have a timescale of < 20 days. When one takes M3 (Fig. 21) 
and M4 (Fig. 22) the situation is also somewhat different. 
The probability of observing a 40-60 days top-hat fluctua- 
tion is now of about 5%, although most features are due to 
a small gap of about 50 days around day 1815 (see Paper I 
and Fig. 2). Finally, Figs. 19-22 inform on the true ability 
of the IAC-80 telescope to detect microlensing fluctuations 
in an observational DLC free from gaps: a "peak" with a 
timescale > 40 days should be interpreted as a feature re- 
lated to microlensing or other mechanisms different to the 
observational noise, while as mainly caused by the poor res- 
olution at the expected amplitudes within the interval [- 50 
mmag, + 50 mmag], the < 20 days microlensing "peaks" 
cannot be resolved. Even in the unlikely case of very short- 
timescale microlensing signal with high amplitude, due to the 
smoothing by both the re-reduction and the binning as well as 
the current uncertainty of one week in the true time delay, it 
would be not possible to reliably reconstruct the microlensing 
"peaks". 



5 CONCLUSIONS 

Several ~ 1 metre class telescopes around the world are at 
present involved in different optical monitoring programs 
of quasars with the goal to detect microlensing. There are 
at least two "modest" telescopes searching for microlensing 
signal related to a far elliptical galaxy (which is responsible, 
in part, for the gravitational mirage Q0957+561). The data 
taken at Whipple Observatory 1.2 m telescope and at Teide 
Observatory IAC-80 telescope together with the photome- 
try from a 3.5 m telescope (at Apache Point Observatory) 
represent a great effort in order to obtain an accurate time 
delay in Q0957+561, follow the long-timescale microlensing 
event in that system and find some evidence in favour of very 
rapid and rapid microlensing (Kundic et al. 1995, 1997; Os- 
coz et al. 1996, 1997; Pijpers 1997; Schild & Thomson 1997; 
Paper I; S96; Pelt et al. 1998; SW98; G98). 

With respect to the very rapid (events with a timescale 
< 3 weeks) and rapid (events with a duration of 1-4 months) 
microlensing, the previous results (before this article) are 
puzzling. The combined photometries (CLCs) from data 
taken at Whipple Observatory only can be well explained 
in the context of a picture including intrinsic variability, 
observational noise and microlensing variability on differ- 
ent timescales: from days to months (e.g., S96). The long- 
timescale microlensing does not play any role in a CLC. In 
particular, S96 reported on the existence of a network of 
rapid events with a few months timescale and an amplitude 
of about ± 50 mmag (these features found by Schild are 
called Schild-events) . However, SW98 concluded that a pic- 
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ture with intrinsic signal and observational noise (without 
any need to introduce very rapid and rapid microlensing) 
is consistent with the observations at Apache Point Obser- 
vatory. SW98 really show a difference light curve in global 
agreement with the zero line, but some doubt remains on the 
ability of the observational noise for producing the negative 
and positive measured events around "peaks" (a "peak" is 
constituted by a set of two or more consecutive deviations 
which have equal sign and are not consistent with zero). In 
any case, SW98 observed no Schild-events. 

In this paper, motivated by the mentioned intriguing re- 
sults on microlensing variability, we analyzed the data from 
our initial monitoring program with the IAC-80 telescope 
(see Paper I). We focused on the possible presence of rapid 
microlensing events in the light curves of QSO 0957+561 and 
the sensitivity of the telescope (using typical observational 
and analysis procedures) to microlensing "peaks". Our con- 
clusions are: 

(i) Using photometric data (in the R band) for the 1996- 
1998 seasons, we made two difference light curves (DLCs). 
The total difference signal, which is based on ~ 1 year of 
overlap between the time-shifted light curve for the A com- 
ponent and the magnitude-shifted light curve for the B com- 
ponent, is in apparent agreement with the absence of mi- 
crolensing signal. We can reject the existence (in our DLCs) 
of events with quarter-year timescale and an amplitude of 
± 50 mmag, and therefore, Schild-events cannot occur al- 
most continuously. On the contrary, they must be either rare 
phenomena (originated by microlensing or another physical 
process) or, because two observatories (Apache Point Ob- 
servatory and Teide Observatory) found no Schild-events, 
untrue fluctuations associated with the observational proce- 
dure and/or the reduction of data at Whipple Observatory. 

(ii) From a very conservative point of view, in our data, 
the amplitude of any hypothetical microlensing signal should 
be in the interval [- 50 mmag, + 50 mmag]. The rms aver- 
ages of the DLCs (global variability) are of about 22 mmag 
(1996/1997 seasons) and 15 mmag (1997/1998 seasons), and 
reasonable constraints on the possible microlensing variabil- 
ity lead to interesting information on the granularity of 
the dark matter in the main lensing galaxy (a cD ellipti- 
cal galaxy) and the size of the source (QSO). Thus the set 
of bounds derived from 1995-1998 seasons (SW98 and this 
work) rules out an important population of MACHOs with 
substellar mass for a small quasar size (Schmidt 1999). 

(iii) In order to settle any doubt on the ability of the ob- 
servational noise for generating the global (rms averages) 
and local (events and other less prominent features) prop- 
erties of the DLCs, we have also carried out several experi- 
ments as "Devil's advocates". The measured variability (the 
rms value, a very rapid event and some minor deviations) in 
the DLC for 1996/1997 seasons can be caused, in a natural 
way, by the observational noise process. In the absence of 
microlensing signal, we proposed two different models (Ml 
and M2; see subsection 3.1) whose associated photometries 
(simulations) are consistent with the observations. In addi- 
tion, the DLC for 1997/1998 seasons is a quasi-featureless 
trend with relatively small rms average. To explain the vari- 
ability in our second observational DLC, we again showed 
that microlensing is not necessary. Two new models (M3 
and M4; see subsection 3.2) only including the reconstruc- 



tion of the intrinsic signal (assumed as a polynomial or a 
correlated stationary process) and a Gaussian observational 
noise process, led to simulated DLCs in agreement with the 
measured behaviour. 

(iv) We finally show that from a typical monitoring with 
our telescope (observing times, method of analysis, etc.) 
is not possible the resolution of < 20 days microlensing 
"peaks". The confusion with noise does not permit the sep- 
aration between true microlensing features and "peaks" due 
to the observational noise. However, all hypothetical "peaks" 
with a timescale > 40 days must be interpreted as phenom- 
ena which are not associated with the observational noise 
(e.g., microlensing fluctuations). At intermediate timescales 
(of about one month) the situation is somewhat intricate. 
Given a measured DLC, the probability of observing one 
noise "peak" (with a duration of about 30 days) is less than 
10%. Therefore, if we search for microlensing signal and find 
an "intermediate peak", the relative probabilities that the 
fluctuation is a noise feature or a microlensing "peak" are 
< 1:10. 
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